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It is known that within the interacting electron model Hamiltonian for the one-dimensional i- 
filled band, the singlet ground state is a Wigner crystal only if the nearest neighbor electron-electron 
repulsion is larger than a critical value. We show that this critical nearest neighbor Coulomb 
interaction is different for each spin subspace, with the critical value decreasing with increasing 
spin. As a consequence, with the lowering of temperature, there can occur a transition from a 
Wigner crystal charge-ordered state to a spin-Peierls state that is a Bond-Charge-Density Wave 
with charge occupancies different from the Wigner crystal. This transition is possible because 
spin excitations from the spin-Peierls state in the ^-filled band are necessarily accompanied by 
changes in site charge densities. We apply our theory to the ^-filled band quasi-one-dimensional 
organic charge-transfer solids in general and to 2:1 tetramethyltetrathiafulvalene (TMTTF) and 
tetramethyltetraselenafulvalene (TMTSF) cationic salts in particular. We believe that many recent 
experiments strongly indicate the Wigner crystal to Bond-Charge-Density Wave transition in several 
members of the TMTTF family. We explain the occurrence of two different antiferromagnetic 
phases but a single spin-Peierls state in the generic phase diagram for the 2:1 cationic solids. The 
antiferromagnetic phases can have either the Wigner crystal or the Bond-Charge-Spin-Density Wave 
charge occupancies. The spin-Peierls state is always a Bond-Charge-Density Wave. 

PACS numbers: 71.30.+h, 71.45.Lr, 74.70.Kn 



I. INTRODUCTION 

Spatial broken symmetries in the quasi-one- 
dimensional (quasi-lD) ^-filled organic charge- 
transfer solids (CTS) have been of strong 
experimentall^^^ao^^^^ and 
theoreticall 8 .^^^^ 2 ^^^ interest. The bro- 
ken symmetry states include charge order (hereafter 
CO, this is usually accompanied by intramolecular 
distortions), intermolecular lattice distortions (hereafter 
bond order wave or BOW), antiferromagnetism (AFM) 
and spin-Peierls (SP) order. Multiple orderings may 
compete or even coexist simultaneously. Interestingly, 
these unconventional insulating states in the CTS are 
often proximate to superconductivity 2 !, the mecha- 
nism of which has remained perplexing after intensive 
investigations over several decades. Unconventional 
behavior at or near ^-filling has also been observed 
in the quasi-two-dimensional organic CTS with higher 
superconducting critical temperature a 11 ! 28 ' 29 ! 30 ' 31 , 
sodium cobaltat o 32 ' 33 i 34 and oxides of titanium 3 ^ and 
vanadiu m 36 i 3 . 

In spite of extensive research on ID instabilities in the 
CTS, detailed comparisons of theory and experiments re- 
main difficult. Strong electron-electron (e-e) Coulomb 
interactions in these systems make calculations particu- 
larly challenging, and with few exceptions 3 ^ 9 - existing 
theoretical discussions of broken symmetries in the in- 
teracting ^-filled band have been limited to the ground 
stat elU&2iL2L2&224tfU&. This leaves a number of im- 
portant questions unresolved, as we point out below. 
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FIG. 1: Schematic of the proposed T-vs.-P phase diagram for 
(TMTCF) 2 X, where C=T or S, along with the charge occu- 
pancies of the sites in the low T phases as determined in this 
work. The P-axis reflects the extent of interchain coupling. 
Open (filled) arrows indicate the ambient pressure locations 
of TMTTF (TMTSF) salts. 



Quasi- ID CTS undergo two distinct phase transitions 
as the temperature (hereafter T) is reduced 40 . The 4kp 
transition at higher T involves charge degrees of free- 
dom (T4fc jr ~ 100 K), with the semiconducting state be- 
low T^k F exhibiting either a dimerized CO or a dimer- 
ized BOW. Charges alternate as 0.5 + e and 0.5 - e 
on the molecules along the stack in the dimerized CO. 
The site charge occupancy in the dimerized CO state is 
commonly written as • • • 1010- • • (with T' and '0' denot- 
ing charge-rich and charge-poor sites), and this state is 
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also referred to as the Wigner crystal. The 4kp BOW 
has alternating intermolecular bond strengths, but site 
charges are uniformly 0.5. It is generally accepted that 
the dimer CO (dimer BOW) is obtained for strong (weak) 
intermolecular Coulomb interactions (the intramolecular 
Coulomb interaction can be large in cither case, see Sec- 
tion III). At T < T 2 k F ~ 10 - 20 K, there occurs a sec- 
ond transition in the CTS, involving spin degrees of free- 
dom, to either an SP or an AFM state. Importantly, 
the SP and the AFM states both can be derived from 
the dimer CO or the dimer BOW. Coexistence of the SP 
state with the Wigner crystal would require that the SP 
state has the structure l = = l- -- 0- -- l, with singlet 
"bonds" alternating in strength between the charge-rich 
sites of the Wigner crystal. Similarly, coexisting AFM 
and Wigner crystal would imply site charge-spin occu- 
pancies | Q J. 0. The occurrence of both have been sug- 
gested in the literatur e 20 ' 26 ' 38 ' 39 . The SP state can also 
be obtained from dimerization of the dimer 4kp BOW, 
in which case there occurs a spontaneous transition from 
the uniform charge-density state to a coexisting bond- 
charge-density wave (BCDW)±&22>25 1 - 1 = • • • 0. The 
unit cells here are dimers with site occupancies '1' and 
'0' or '0' and '1', with strong and weak interunit 1-1 and 
0- • • bonds, respectively. The AFM state with the same 
site occupancies • ■ • 1100- • • is referred to as the bond- 
charge-spin density wave (BCSDW) 21 and is denoted as 

Tl oo. 

The above characterizations of the different states and 
the related bond, charge and spin patterns are largely 
based on ground state calculations. The mechanism of 
the transition from the 4kp state to the 2kp state lies 
outside the scope of such theories. Questions that remain 
unresolved as a consequence include: (i) Does the nature 
of the 4kp state (charge versus bond dimerized) predeter- 
mine the site charge occupancies of the 2kp state? This 
is assumed in many of the published works, (ii) What de- 
termines the nature of the 2kp state (AFM versus SP)? 
(iii) How does one understand the occurrence of two dif- 
ferent AFM phases (hereafter AFM1 and AFM2) strad- 
dling a single SP phase in the proposed T vs. P (where P 
is pressure) phase diagram (see Fig. [T]^ 9 - for the cationic i- 
filled band quasi-lD CTS? (iv) What is the nature of the 
intermediate T states between the 4kj? dimerized state 
and the 2kp tetramerized state? As we point out in the 
next section, where we present a brief review of the ex- 
perimental results, answers to these questions are crucial 
for a deeper understanding of the underlying physics of 
the i-filled band CTS. 

In the present paper we report the results of our calcu- 
lations of T-dependent behavior within theoretical mod- 
els incorporating both e-e and electron-phonon (e-p) in- 
teractions to answer precisely the above questions. One 
key result of our work is as follows: in between the strong 
and weak intermolecular Coulomb interaction parame- 
ter regimes, there exists a third parameter regime within 
which the intermolecular Coulomb interactions are in- 
termediate, and within which there can occur a novel 



transition from a Wigner crystal CO state to a BCDW 
SP state as the T is lowered. Thus the charge or bond 
ordering in the 4k F phase does not necessarily decide 
the same in the 2kp state. For realistic intramolecular 
Coulomb interactions (Hubbard U), we show that the 
width of this intermediate parameter regime is compara- 
ble to the strong and weak intersite interaction regimes. 
We believe that our results are directly applicable to the 
family of the cationic CTS (TMTTF) 2 X, where a redis- 
tribution of charge upon entering the SP state from the 
CO state has been observecU^i 6 - in X=AsF 6 and PF 6 . A 
natural explanation of this redistribution emerges within 
our theory. The SP state within our theory is unique and 
has the BCDW charge occupancy, while the two AFM re- 
gions in the phase diagram of Fig. [1] have different site 
occupancies. Our theory therefore provides a simple di- 
agnostic to determine the pattern of CO coexisting with 
low-temperature magnetic states in the ^-filled CTS. 

In addition to the above theoretical results directly 
pertaining to the quasi- ID CTS, our work gives new in- 
sight to excitations from a SP ground state in a non- 
half-filled band. In the case of the usual SP transition 
within the 1-filled band, the SP state is bond-dimerized 
at T = and has uniform bonds at T > T%k F ■ The site 
charges are uniform at all T. This is in contrast to the 
1-filled band, where the SP state at T = is bond and 
charge-tetramerized and the T > T 2 k F state is dimerized 
as opposed to being uniform. Furthermore, we show that 
the high T phase here can be either charge- or bond- 
dimerized, starting from the same low T state. This 
clearly requires two different kinds of spin excitations in 
the ^-filled band. We demonstrate that spin excitations 
from the SP state in the 1-filled band can lead to two 
different kinds of defects in the background BCDW. 

In the next section we present a brief yet detailed sum- 
mary of relevant experimental results in the quasi- ID 
CTS. The scope of this summary makes the need for hav- 
ing T-dependent theory clear. Following this, in Section 
III we present our theoretical model along with conjec- 
tures based on physical intuitive pictures. In Section IV 
we substantiate these conjectures with accurate quantum 
Monte Carlo (QMC) and exact diagonalization (ED) nu- 
merical calculations. Finally in Section V we compare 
our theoretical results and experiments, and present our 
conclusions. 



II. REVIEW OF EXPERIMENTAL RESULTS 

Examples of both CO and BOW broken symmetry at 
T < T 4kF are found in the ^-filled CTS. The 4k F phase 
in the anionic 1:2 CTS is commonly bond-dimerized. The 
most well known example is MEM(TCNQ)2, which un- 
dergoes a metal-insulator transition accompanied with 
bond-dimerization at 335 K— . Site charges are uniform 
in this 4kp phase. The 2kp phase in the TCNQ-based 
systems is universally SP and not AFM. The SP transi- 
tion in MEM(TCNQ) 2 occurs below T 2kF = 19 K, and 
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low T neutron diffraction measurements of deuterated 
samples- 1 have established that the bond tetramerization 
is accompanied by 2kp CO ■ • ■ f fOO- ■ • . X-ra y 42 ' 43 and 
neutron diffraction 44 experiments have confirmed a sim- 
ilar low T phase in TEA(TCNQ) 2 . We will not discuss 
these further in the present paper, as they are well de- 
scribed within our previous wor k 18 ' 25 . We will, however, 
argue that the SP ground state in (DMe-DCNQl) 2 Ag 4 ^ 
(as opposed to AFM) indicates • ■ • 1 fOO- ■ • CO in this. 

The cationic (TMTCF) 2 X, C= S and Se, exhibit more 
variety, presumably because the counterions affect pack- 
ing as well as site energies in the cation stack. Differ- 
ences between systems with centrosymmetric and non- 
centrosymmetric anions are also observed. Their overall 
behavior is summarized in Fig. [1] where as is custom- 
ary pressure P can also imply larger interchain coupling. 
We have indicated schematically the possible locations of 
different materials on the phase diagram. The most sig- 
nificant aspect of the phase diagram is the occurrence of 
two distinct antiferromagnetic phases, AFMf and AFM2, 
straddling a single SP phase. 

Most TMTTF lie near the low P region of the phase 
diagram and are insulating already at or near room tem- 
perature because of charge localization, which is due to 
the intrinsic dimerization along the cationic stack o 1 ' 4 ' 11 . 
CO at intermediate temperatures Tco has been found 
in dielectric permittivity 4 , XMlii and ESB4 6 - experi- 
ments on materials near the low and intermediate P 
end. Although the pattern of the CO has not been de- 
termined directly, the observation of ferroelectric behav- 
ior below Tco is consistent with • • • f Of 0- • • type CO 
in this region 5 -^. With further lowering of T, most 
(TMTTF) 2 X undergo transitions to the AFMf or SP 
phase (with X = Br a possible exception, see below). 
X = SbFg at low T lies in the AFMf region^, with a very 
high Tco and relatively low Neel temperature Tjv = 8 
K. As the schematic phase diagram indicates, pressure 
suppresses both Tco and Tjv in this region. For P > 0.5 
GPa, (TMTTF) 2 SbF 6 undergoes a transition from the 
AFMf to the SP phased, the details of which are not com- 
pletely understood; any charge disproportionation in the 
SP phase is smal£. (TMTTF) 2 Re0 4 also has a relatively 
high Tco — 225 K, but the low T phase here, reached 
following an anion-ordering transition is spin singlelji 4 -. 
Nakamura et al. have suggested, based on NMR exper- 
iments, that the CO involves the Wigner crystal state, 
but the low T state is the • • • ilOO- • • BCDWil Further 
along the P axis lie X = AsF 6 and PF 6 , where Tco are re- 
duced to f 00 K and 65 K, respectively^. The low T phase 
in both cases is now SP. Neutron scattering experiments 
on (TMTTF) 2 PF 6 have found that the lattice distortion 
in the SP state is the expected 2kp BOW distortion, 
but that the amplitude of the lattice distortion is much 
smaller— than that found in other organic SP materials 
such as MEM(TCNQ) 2 . The exact pattern of the BOW 
has not been determined yet. Experimental evidence ex- 
ists that some form of CO persists in the magnetic phases. 
For example, the splitting in vibronic modes below Tco 



in (TMTTF) 2 PF 6 and (TMTTF) 2 AsF 6 , a signature of 
charge disproportionation, persists into the SP phase 17 , 
indicating coexistence of CO and SP. At the same time, 
the high T CO is in competition with the SP ground 
stated, as is inferred from the different effects of pressure 
on Tco and T$p'- while pressure reduces Ted h in- 
creases Tsp- This is in clear contrast to the effect of pres- 
sure on Tjv in X = SbFg . Similarly, deuteration of the hy- 
drogen atoms of TMTTF increases Tco but decreases 6 ^ 
Tsp- That higher Tco is accompanied by lower T$p for 
centrosymmetric X (T SP = 16.4 K in X = PF 6 and If .1 K 
in X = AsFg) has also been noted 4 ^. This trend is in ob- 
vious agreement with the occurrence of AFM instead of 
SP state under ambient pressure in X = SbF 6 . Most in- 
terestingly, Nakamura et al. have very recently observed 
redistribution of the charges on the TMTTF molecules 
in (TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 as these systems 
enter the SP phase from CO state o 15 ' 16 . Charge dispro- 
portionation, if any, in the SP phase is much smaller 
than in the CO phas o 15 i 16 , which is in apparent agree- 
ment with the above observations^ 4 in X = Re04 and 
SbF 6 . 

The bulk of the (TMTTF) 2 X therefore lie in the 
AFMf and SP regions of Fig. [Q (TMTSF) 2 X, in 
contrast, occupy the AFM2 region. Coexisting 2kp 
CDW and spin-density wave, SDW, with the same 
2kp periodicit y 49 i 50 here is explained naturally as the 
• • • 1100- • • BCSDWSii 2 ^ 1 .. fn contrast to the TMTTF 
salts discussed above, charge and magnetic ordering in 
(TMTTF) 2 Br occur almost simultaneousl y 46 ' 52 . X-ray 
studies of lattice distortions point to similarities with 
(TMTSF) 2 PF 9 4 ^, indicating that (TMTTF) 2 Br is also 
a • • • If 00- • • BCSDWSi. We do not discuss AFM2 re- 
gion in the present paper, as this can be found in our 
earlier wor k 21 ' 25 . 



III. THEORETICAL MODEL AND 
CONJECTURES 

The ID Hamiltonian we investigate is written as 

H = Hssh + HhoI + H ee (la) 
Hssh = t^[l + a(a\ + ai)](cl a c l+1 . a + h.c.) 

i 

+ hhj s ^^a\ai (lb) 

i 

H H oi = g^ibl + b^m + HujH^blbi (lc) 

i i 

n i,t n i,l + V /; n i n i+l (l d ) 
i i 

In the above, cj creates an electron with spin a (T,|) on 
molecular site i, 7ij )Cr = c\ a .c^ rJ is the number of electrons 
with spin a on site i, and rij = Y] n^o-- U and V are the 
on-site and intersite Coulomb repulsions, and a\ and b\ 
create (dispersionless) Su-Schrieffer-Heeger (SSH) 5 - 3 - and 
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Holstein (Hoi) 54 phonons on the ith bond and site respec- 
tively, with frequencies u>$ and u>h- Because the Peierls 
instability involves only phonon modes near q = n, keep- 
ing single dispersionless phonon modes is sufficient for 
the Peierls transitions to occur— i 5 - 6 -. Although purely ID 
calculations cannot yield a finite temperature phase tran- 
sition, as in all low dimensional theorie s 57 ! 58 we antici- 
pate that the 3D ordering in the real system is principally 
determined by the dominant ID instability. 

The above Hamiltonian includes the most important 
terms necessary to describe the family of quasi-lD CTS, 
but ignores nonessential terms that may be necessary for 
understanding the detailed behavior of individual sys- 
tems. Such nonessential terms include (i) the intrinsic 
dimcrization that characterizes many (TMTTF^X, (ii) 
interaction between counterions and the carriers on the 
quasi-lD cations stacks, (iii) interchain Coulomb inter- 
action, and (iv) interchain hopping. Inclusion of the in- 
trinsic dimerization will make the Wigner crystal ground 
state even less likely 2 ^, and this is the reason for exclud- 
ing it. We have verified the conclusions of reference 
from exact diagonalization calculations. The inclusion of 
interactions with counterions may enhance the Wigner 
crystal ordering 5 - 2 ^, for some (TMTTF) 2 X. We will dis- 
cuss this point further below, and argue that it is im- 
portant in the AFM1 region of the phase diagram. The 
effects of intrinsic dimerization and counterion interac- 
tions can be reproduced by modifying the V/\t\ in our 
Hamiltonian, and thus these are not included explicitly. 
Rather the V in Eq. (|la|) should be considered as the 
effective V for the quasi-lD CTS. Interchain hopping, at 
least in the TMTTF (though not in the TMTSF), are in- 
deed negligible. The interchain Coulomb interaction can 
promote the Wigner crystal within a rectangular lattice 
framework, but for the realistic nearly triangular lattice 
will cause frustration, thereby making the BCDW state 
of interest here even more likely. We therefore believe 
that Hamiltonian captures the essential physics of 
the quasi-lD CTS. 

For applications to the quasi-lD CTS we will be inter- 
ested in the parameter regim o 18 i 24 i 25 ' 26 \t\ = 0.1 — 0.25 
eV, U/\t\ — 6 — 8. The exact value of V is less known, 
but since the same cationic molecules of interest as well 
as other related molecules (for e.g., HMTTF, HMTSF, 
etc.) also form quasi-lD ^-filled band Mott-Hubbard 
semiconductors with short-range antiferromagnetic spin 
correlation a 59 ! 60 , it must be true that V < \U (since for 
V > \U the ID i-filled band is a CDW 6 ^). Two other 
known theoretical results now fix the range of V that 
will be of interest. First, the ground state of the ^-filled 
band within Hamiltonian (|la[) in the limit of zero e-p 
coupling is the Wigner crystal ■ ■ • 1010- • ■ only for suffi- 
ciently large V > V C (U). Second, V C {U ->■ oo) = 2\t\, and 
is larger for finite £ / 19 i 24 i 25 . With the known U/\t\ and 
the above restrictions in place, it is now easily concluded 
that (i) the Wigner crystal is obtained in the CTS for a 
relatively narrow range of realistic parameters, and (ii) 
even in such cases the material's V is barely larger than 
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We now go beyond the above ground state theories 
of spatial broken symmetries to make the following cru- 
cial observation: each different spin subspace of Hamil- 
tonian Ua]) must have its own V c at which the Wigner 
crystal is formed. This conclusion of a spin- dependent 
V c = V C (U,S) follows from the comparison of the fer- 
romagnetic subspace with total spin S = S max and the 
S = subspace. The ferromagnetic subspace is equiva- 
lent to the ^-filled spinless fermion band, and therefore 
V C (U, S max ) is independent of U and exactly 2\t\. The 
increase of V C (U, S = 0) with decreasing i s then 

clearly related to the occurrence of doubly occupied and 
vacant sites at finite U in the S — subspace. Since 
the probability of double occupancy (for fixed U and V) 
decreases monotonically with increasing 5, we can in- 
terpolate between the two extreme cases of S — and 
S = S max to conclude V C (U, S) > V C (U,S+1). We prove 
this explicitly from numerical calculations in the next 
section. 

Our conjecture regarding spin-dependent V C (U) in turn 
implies that there exist three distinct parameter regimes 
for realistic U and V: (i) V < V{S max ) = 2\t\, in 
which case the ground state is the BCDW with 2kp CO 
and the high temperature 4k^ state is a BO W 18 ' 25 ; (ii) 
V > V C (U, S = 0), in which case both the 4kp and 
the 2kp phases have Wigner crystal CO; (iii) and the 
intermediate regime 2\t\ < V < V C (U, S = 0). Pa- 
rameter regime (i) has been discussed in our previous 
wor k 18 i 21 ' 25 i 63 . This is the case where the description of 
the ^-filled band as an effective ^-filled band is appropri- 
ate: the unit cell in the 4kp state is a dimer of two sites, 
and the 2kp transition can be described as the usual SP 
dimerization of the dimer lattice. We will not discuss 
this parameter region further. Given the U/\t\ values for 
the CTS, parameter regime (ii) can have rather narrow 
width. For U/\t\ = 6, V C (U,S = 0) = 3.3|i|, and there 
is no value of realistic V for which the ground state is a 
Wigner crystal. For U/\t\ = 8, V C (U, S = 0) = 3.0|t|, and 
the widths of parameter regime (iii), 2\t\ < V < 3\t\ and 
of parameter regime (ii), 3|£| < V < 4\t\, are comparable. 

We will discuss parameter regime (ii) only briefly, as 
the physics of this region has been discussed by other 
author s 20 ! 26 . We investigated thoroughly the intermedi- 
ate parameter regime (iii), which has not been studied 
previously. Within the intermediate parameter regime 
we expect a novel transition from the BCDW in the 2k^? 
state at low T to a • • ■ 1010- • • CO in the 4kp phase at 
high T that has not been discussed in the theoretical lit- 
erature. 

Our observation regarding the T-dependent behavior 
of these CTS follows from the more general observation 
that thermodynamic behavior depends on the free en- 
ergy and the partition function. For the strong e-e inter- 
actions of interest here, thermodynamics of ID systems 
at temperatures of interest is determined almost entirely 
by spin excitations. Since multiplicities of spin states 
increase with the total spin S, the partition function is 
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dominated by high (low) spin states with large (small) 
multiplicities at high (low) temperatures. While at T=0 
such a system must be a BCDW for V < V C (U, S = 0), as 
the temperature is raised, higher and higher spin states 
begin to dominate the free energy, until V for the ma- 
terial in question exceeds V C (U,S), at which point the 
charge occupancy reverts to • • • 1010- • • We demonstrate 
this explicitly in the next section. A charge redistribu- 
tion is expected in such a system at T2k F , as the devi- 
ation from the average charge of 0.5 is much smaller in 
the BCDW than in the Wigner crystal^. 

The above conjecture leads to yet another novel impli- 
cation. The ground state for both weak and intermediate 
intersite Coulomb interaction parameters are the BCDW, 
even as the 4kp phases are different in the two cases: 
BOW in the former and CO in the latter. This necessarily 
requires the existence of two different kinds of spin exci- 
tations from the BCDW. Recall that within the standard 
theory of the SP transition in the |-filled band^^^, 
thermal excitations from the T = ground state gener- 
ates spin excitations with bond-alternation domain walls 
(solitons), with the phase of bond alternation in between 
the solitons being opposite to that outside the solitons. 
Progressive increase in T generates more and more soli- 
tons with reversed bond alternation phases, until overlaps 
between the two phases of bond alternations lead ulti- 
mately to the uniform state. A key difference between 
the i-filled and ^-filled SP states is that site charge oc- 
cupancies in spin excitations in the former continue to be 
uniform, while they are necessarily nonuniform in the lat- 
ter case, as we show in the next section. We demonstrate 
that defect centers with two distinct charge occupancies 
(that we will term as type I and II) , depending upon the 
actual value of V, are possible in the ^-filled band. Pre- 
ponderance of one or another type of defects generates 
the distinct 4k_F states. 



cal ensemble. Quantum phonons are treated within SSE 
by directly adding the bosonic phonon creation and an- 
nihilation operators in Eq. (fTaj) to the series expansion 67 . 
An upper limit in the phonon spectrum must be imposed, 
but can be set arbitrarily large to avoid any systematic 
errors^. For the results shown below we used a cutoff of 
100 SSH phonons per bond and either 30 (for g = 0.5) or 
50 (for g — 0.75) Holstein phonons per site. 

The observables we calculate within SSE are the 
standard wave-vector dependent charge structure factor 
S p (q), defined as, 

j.k 

and charge and bond-order susceptibilities Xp{l) an d 
Xb(q), defined as, 

XM = ij2 ei9ij ' k) f 1 dr(O](r)OU0)) (3) 

In Eqs. ([2|) and (J3j> N is the number of lattice sites, Oj = 

n i,T + °f = E<t( c ]+i^ c .j> + h.c), and (3 is the 
inverse temperature in units of t. 

The presence of CO or BOW can be detected by the 
divergence of the 2\ip or 4k p charge or bond-order sus- 
ceptibility as a function of increasing system size. Strictly 
speaking, in a purely ID model these functions diverge 
only at T = 0; as already explained above, we make the 
reasonable assumption— that in the presence of realis- 
tic inter-chain couplings transitions involving charge or 
bond-order instabilities, as determined by the dominant 
susceptibility, occur at finite T. 

A. Spin-dependent V C (U) 



IV. RESULTS 

We present in this section the results of QMC investi- 
gations of Eq. (|la|). for both zero and nonzero e-p cou- 
plings, and of ED studies of the adiabatic (semiclassical) 
limit of Eq. ([Taj) . Using QMC techniques, we demon- 
strate explicitly the spin-dependence of V C (U), as well as 
the transition from the Wigner crystal to the BCDW for 
the intermediate parameter regime (iii). The ED studies 
demonstrate the exotic nature of spin excitations from 
the BCDW ground state. In what follows, all quantities 
are expressed in units of \t\ (|t|=l). 

The QMC method we use is the Stochastic Series Ex- 
pansion (SSE) method using the directed loop update 
for the electron degrees of freedom 6 - 6 . For ID fermions 
with nearest-neighbor hopping SSE provides statistically 
exact results with no systematic errors. While the SSE 
directed-loop method is grand canonical (with fluctuat- 
ing particle density), we restrict measurements to only 
the i-filled density sector to obtain results in the canoni- 



We first present computational results within Hamil- 
tonian (|la[) in the absence of e-p coupling to demon- 
strate that the V c (U) at which the Wigner crystal order 
is established in the lowest state of a given spin sub- 
space decreases with increasing spin S. Our computa- 
tional approach conserves the total z-component of spin 
S z and not the total S. Since the Lieb-Mattis theorem 6 ^ 
E(S) < E(S+1), where E(S) is the energy of the lowest 
state in the spin subspace S, applies to the ID Hamilto- 
nian (| 1 a|) . and since in the absence of a magnetic field 
all S z states for a given S are degenerate, our results for 
the lowest state within each different S z must pertain to 
S = S Z . 

To determine V C (U, S) we use the fact that the purely 
electronic model is a Luttinger liquid (LL) for V < V c 
with correlation functions determined by a single expo- 
nent K p (see Reference HH for a review). The Wigner 
crystal state is reached when K p — \. The exponent 
Kp may be calculated from the long-wavelength limit of 
S p (q) 70 . In Fig. j2ja) we have plotted our calculated K p 
for U — 8 as a function of V for different S z sectors. The 
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FIG. 2: V c for Eq. (|la[) in the limit of zero e-p interactions 
(a = g — 0) as a function of &. Results are for a iV=32 
site periodic ring with U — 8. For iV = 32 S z ~8 corresponds 
to the fully polarized (spinless fermion) limit, (a) Luttinger 
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temperature chosen is small enough ((3 = 27V) that in all 
cases the results correspond to the lowest state within 
a given S z . In Fig. [2{b) we have plotted our calculated 
V C (U — 8), as obtained from Fig. HJa), as a function of 
S z . V c is largest for S z — and decreases with increas- 
ing S z , in agreement with the conjecture of Section III. 
Importantly, the calculated V c for S z — 8 is close to the 
correct limiting value of 2, indicating the validity of our 
approach. We have not performed any finite-size scaling 
in Fig. [2j which accounts for the the slight deviation from 
the exact value of 2. 
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FIG. 3: (color online) QMC results for the temperature- 
dependent charge susceptibilities for a N=64 site periodic 
ring with U = 8, V = 2.75, a = 0.15, lo s = 0.1, g = 0.5, 
and ujh = 0.5. (a) and (b) Wavevector-dependent charge 
and bond-order susceptibilities, (c) 2ki? and 4ki? charge sus- 
ceptibilities as a function of temperature, (d) 2kF and 4kF 
bond-order susceptibilities as a function of temperature. If 
error bars are not shown, statistical error bars are smaller 
than the symbol sizes. Lines are guides to the eye. 



likelihood of the Wigner crystal CO. We report results 
only for intermediate and strong weak intersite Coulomb 
interactions; the weak interaction regime V < 2 has been 
discussed extensively in our previous work2£. 



1. 2 < V < V C (U,S = 0) 



B. T-dependent susceptibilities 

We next present the results of QMC calculations within 
the full Hamiltonian (fl~a|) . To reproduce correct relative 
energy scales of intra- and intermolecular phonon modes 
in the CTS, we choose luh > wj, specifically loh = 0.5 
and u)s = 0.1 in our calculations. Small deviations from 
these values do not make any significant difference. In 
all cases we have chosen the electron-molecular vibration 
coupling g larger than the coupling between electrons and 
the SSH phonons a, thereby deliberately enhancing the 



Our results are summarized in Figs. [3] [5J where we re- 
port results for two different V of intermediate strength 
and several different e-p couplings. Fig. [3] first shows re- 
sults for relatively weak SSH coupling a. The charge sus- 
ceptibility is dominated by a peak at 4ki? = ir (Fig.[3ja)). 
Figs. 03c) and (d) show the T-dependence of the 2kp 
as well as 4k^? charge and bond susceptibility. For 
V < V C (U,S = 0) the ikp charge susceptibility does 
not diverge with system size, and the purely ID system 
remains a LL with no long range CO at zero tempera- 
ture. The dominance of Xp(4&f) over Xp(^f) suggests, 
however, that • • ■ 1010- ■ • CO will likely occur in the 3D 
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FIG. 4: (color online) Same as Fig. [3] but with parameters 
U = 8, V = 2.25, a = 0.27, lu s = 0.1, and lo h = 0.5. In panels 
(a) and (b), data are for g — 0.50 only. In panels (c) and (d), 
data for both g — 0.50 and g = 0.75 are shown. Arrows 
indicate temperature where Xp(2&f) = Xp(4&f) (solid and 
broken arrows correspond to g — 0.50 and 0.75, respectively.) 



system, especially if this order is further enhanced due to 
interactions with the counterions 9 -. We have plotted the 
bond susceptibilities xsi^kp) and xb(4&f) in Fig.^d). 
A SP transition requires that XB(2kp) diverges as T — > 0. 
XB^kp) weaker than XB(4fcp ) at low T indicates that 
the SP order is absent in the present case with weak SSH 
e-p coupling. This result is in agreement with earlier 
result^ that the SP order is obtained only above a critical 
a c (a c may be smaller for the infinite system than found 
in our calculation). The most likely scenario with the 
present parameters is the persistence of the ■ • • 1010- • • 
CO to the lowest T with no SP transition. These pa- 
rameters, along with counterion interactions, could then 
describe the (TMTTF) 2 X materials with an AFM ground 
state 9 . 

In Figs. Q] and Fig. we show our results for larger 
SSH c-p coupling a and two different intersite Coulomb 
interaction V—2.25 and 2.75. The calculations of Fig. Q] 
were done for two different Holstein e-p couplings g. For 
both the V parameters, x p (4fcp) dominates at high T 
but x p (2fci?) is stronger at low T in both Fig. 0Jc) and 
Fig. \5[c) . The crossing between the two susceptibilities 
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FIG. 5: (color online) Same as Fig. [3] but with parameters 
U = 8, V = 2.75, a = 0.27, uj s = 0.1, g = 0.5, and uj h = 0.5. 
Arrow indicates temperature where Xp(2&f) = Xp(4&f). 



is clear indication that in the intermediate parameter 
regime, as T is lowered the 2kp CDW instability domi- 
nates over the 4kp CO. 

The rise in Xp(2fcj?) at low T is accompanied by a 
steep rise in xsi^kp) in both cases (see Fig. @|d) and 
Fig.[5Jd)). Importantly, unlike in Fig. [3^d), xsi^kp) in 
these cases clearly dominates over xs(4/cp) by an order 
of magnitude at low temperatures. There is thus a clear 
signature of the SP instability for these parameters. The 
simultaneous rise in XB^kp) and x P (2/cp) indicates that 
the SP state is the • • ■ 1100- • ■ BCDW. Comparison of 
Figs. 2] and [5] indicates that the effect of larger V is to 
decrease the T where the 2kp and 4kp susceptibilities 
cross. Since larger V would imply larger Tco, this result 
implies that larger Tco is accompanied by lower Tsp- 
Our calculations are for relatively modest a < g. Larger 
a (not shown) further strengthens the divergence of 2k p 
susceptibilities. 

The motivation for performing the calculations of 
Fig. 4 with multiple Holstein couplings was to deter- 
mine whether it is possible to have a Wigner crystal at 
low T even for V < V(U,S — 0), by simply increas- 
ing g. The argument for this would be that in strong 
coupling, increasing g amounts to an effective increase 
in V— . The Holstein coupling cannot be increased arbi- 
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FIG. 6: (color online) Same as Fig. [3j but with parameters 
U = 8, V = 3.5, a = 0.24, lo s = 0.1, g = 0.5, and lu h = 0.5. 



trarily, however, as beyond a certain point, g promotes 
formation of on-site bipolarons^. Importantly, the co- 
operative interaction between the 2kp BOW and CDW 
in the BCDW2i>25 implies that in the V < V(U, 5 = 0) 
region, larger g not only promotes the 4kp CO but also 
enhances the BCDW. In our calculations in Fig. H£b) 
and (c) we we find both these effects: a weak increase in 
Xp(4/cf) at intermediate T, and an even stronger increase 
in the T — > values of Xp(2&f) and xs(2kp). Actually, 
this result is in qualitative agreement with our obser- 
vation for the ground state within the adiabatic limit 
of Eq. (la} that in the range < V < V(U,S = 0), 
V enhances the BCDW—. The temperature at which 
X P (2kp) and Xp(4&f) cross does not change significantly 
with larger g. Our results for g = 0.75 for V — 2.75 (not 
shown) are similar, except that the data are more noisy 
now (probably because all parameters in Fig. 5 are much 
too large for the larger g). We conclude therefore that in 
the intermediate V region, merely increasing g does not 
change the BCDW nature of the SP state. For g to have 
a qualitatively different effect, V should be much closer 
to V(U, S = 0) or perhaps even larger. 



2. V>V c (U,S = 0) 

In principle, calculations of low temperature instabili- 
ties here should be as straightforward as the weak inter- 
site Coulomb interaction V < 2t regime. The 4kp CO - 
AFM1 state t | would occur naturally for the case 
of weak e-p coupling. Obtaining the 4kp CO-SP state 
1 = = 1 • • • ■ • • 1, with realistic V < |Z7 is, however, 
difficult^. Previous work-22, for example, finds this state 
for V > t^U. Recent T-dependent mean-field calculations 
of Seo et aZ.— also fail to find this state for nonzero V. 
There are two reasons for this. First, the spin exchange 
here involves charge-rich sites that are second neighbors, 
and are hence necessarily small. Energy gained upon 
alternation of this weak exchange interactions is even 
smaller, and thus the tendency to this particular form 
of the SP transition is weak to begin with. Second, this 
region of the parameter space involves either large U (for 
e.g., U = 10, for which V C (U,S = 0) ~ 2) or relatively 
large V (for e.g. V C (U, S = 0) ~ 3 for U = 8). In either 
case such strong Coulomb interactions make the applica- 
bility of mean-field theories questionable. 

Our QMC calculations do find the tendency to SP in- 
stability in this parameter region. In Fig. [S] we show 
QMC results for V > V C (U = 8,5 = 0). In contrast to 
Figs. [4] and O x p (4/cf) now dominates over x p (2fci?) at 
all T. The weaker peak at q = 2k p at low T is due to 
small differences in site charge populations between the 
charge-poor sites of the Wigner crystal that arises upon 
bond distortion^., and that adds a small period 4 compo- 
nent to the charge modulation. The bond susceptibility 
Xb [q) has a strong peak at 2kp and a weaker peak at 4kp, 
exactly as expected for the 4kp CO-SP state. Previous 
work has shown that the difference in charge densities 
between the charge-rich and charge-poor sites in the 4k^ 
CO-SP ground state is considerably larger than in the 
BCDW^. 



C. Spin excitations from the BCDW 

The above susceptibility calculations indicate that in 
the intermediate V region (Figs. |4] and [5]) the BCDW 
ground state with • • • 1100- • • CO can evolve into the 
• • • 1010- • • CO as T increases. Our earlier work had 
shown that for weak V, the BCDW evolves into the 4k^ 
BOW at high Ti&2£. Within the standard theory of the 
SP transition££i£ii££, applicable to the i-filled band, the 
SP state evolves into the undistorted state for T > Tsp- 
Thus the 4k p distortions, CO and BOW, take the role of 
the undistorted state in the ^-filled band, and it appears 
paradoxical that the same ground state can evolve into 
two different high T states. We show here that this is inti- 
mately related to the nature of the spin excitations in the 
i-filled BCDW. Spin excitations from the conventional 
SP state leave the site charge occupancies unchanged. 
We show below that not only do spin excitations from 
the i-filled BCDW ground state lead to changes in the 



9 



site occupancies, two different kinds of site occupancies 
are possible in the localized defect states that characterize 
spin excited states here. We will refer to these as type 
I and type II defects, and depending on which kind of 
defect dominates at high T (which in turn depends on 
the relative magnitudes of V and the e-p couplings), the 
4kp state is either the CO or the BOW. 

We will demonstrate the occurrence of type I and II 
defects in spin excitations numerically. Below we present 
a physical intuitive explanation of this highly unusual 
behavior, based on a configuration space picture. Very 
similar configuration space arguments can be found else- 
where in our discussion of charge excitations from the 
BCDW in the interacting ^-filled band^. 

We begin our discussion with the standard ^-filled 
band, for which the mechanism of the SP transition is 
well understoo d 5 ^ 64 ' 65 . Fig. [7(a) shows in valence bond 
(VB) representation the generation of a spin triplet from 
the standard ^-filled band. Since the two phases of bond 
alternation are isoenergetic, the two free spins can sepa- 
rate, and the true wavefunction is dominated by VB dia- 
grams as in Fig. [7(b) , where the phase of the bond alter- 
nation in between the two unpaired spins (spin solitons) 
is opposite to that in the ground state. With increasing 
temperature and increasing number of spin excitations 
there occur many such regions with reversed bond al- 
ternations, and overlaps between regions with different 
phases of bond alternations leads ultimately to the uni- 
form state. 

The above picture needs to be modified for the dimer- 
ized dimer BCDW state in the ^-filled band, in which 
the single site of the i-filled band system is replaced by 
a dimer unit with site populations 1 and (or and 1), 
and the stronger interdimer 1-1 bond (weaker interdimer 
0- ■ • bond) corresponds to the strong (weak) bond in the 
standard SP case. Fig. [7(c), showing triplet generation 
from the BCDW, is the ^-filled analog of Fig. [7(a): a 
singlet bond between the dimer units has been broken to 
generate a localized triplet. The effective repulsion be- 
tween the free spins of Fig. [7(a) is due to the absence 
of binding between them, and the same is expected in 
Fig. [71c). Because the site occupancies within the dimer 
units are nonuniform now, the repulsion between the 
spins is reduced from changes in intradimer as well as 
interdimer site occupancies: the site populations within 
the neighboring units can become uniform (0.5 each), or 
the site occupancies can revert to 1001 from 0110. There 
is no equivalent of this step in the standard SP case. 
The next steps in the separation of the resultant defects 
are identical to that in Fig. [7(b), and we have shown 
these possible final states in Fig. [7(d) and Fig. [7(e), for 
the two different intraunit spin defect populations. For 
V < V C (U, S), defect units with site occupancies 0.5 occu- 
pancies (type I defects) are expected to dominate; while 
for V > V C (U, S) site populations of 10 and 01 (type II 
defects) dominate. From the qualitative discussions it is 
difficult to predict whether the defects are free, in which 
case they are solitons, or if they are bound, in which case 
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FIG. 7: (a) and (b) S = 1 VB diagrams in the Heisenberg SP 
chain. The SP bond order in between the S = | solitons in (b) 
is opposite to that elsewhere, (c) ^-filled band equivalent of 
(a); the singlet bonds in the background BCDW are between 
units containing a pair of sites but a single electron (see text), 
(d) and (e) ^-filled band equivalents of (b). Defect units with 
two different charge distributions are possible. The charge 
distribution will depend on V. 



two of them constitute a triplet. What is more relevant 
is that type I defects generate bond dimerization locally 
(recall that the 4kp BOW has uniform site charge densi- 
ties) while type II defects generate local site populations 
1010, which will have the tendency to generate the 4k^? 
CO. 

The process by which the BCDW is reached from the 
4k^ BOW or the CO as T is decreased from T 2 k F is ex- 
actly opposite to the discussion of spin excitations from 
the BCDW in Fig. [7j Now •••1100--- domain walls 
appear in the background bond-dimerized or charge- 
dimerized state. In the case of the ■ • ■ 1010- • • 4kp state, 
the driving force for the creation of such a domain wall 
is the energy gained upon singlet formation (for V < 
V C (U,S)). 

Fig. [7(e) suggests the appearance of localized 
• • • 1010- • • regions in excitations from the • • • 1100- • • 
BCDW SP state for 2\t\ < V < V C (U,S = 0). We have 
verified this with exact spin-dependent calculations for N 
= 16 and 20 periodic rings with adiabatic e-p couplings, 



H SSH = t^i 1 ~ a ( u i ~ u i+l)}(4,*' 

i 

i 

H Hoi =9^2 v ini + ]-K H vf 



i+l,<7 



h.C.) 



(4a) 
(4b) 



and H ee as in Eq. (lldjl . In the above Ui is the displace- 
ment from equilibrium of a molecular unit and Vi is a 
molecular mode. Note that due to the small sizes of 
the systems we are able to solve exactly, the e-p cou- 
plings in Eq. (|4aj) and Eq. (|4bj) cannot be directly com- 
pared to those in Eq. (TTa[) . In the present case, we 
take the smallest e-p couplings necessary to generate the 
BCDW ground state self-consistently within the adia- 
batic Hamiltonia n 22 i 23 i 25 . We then determine the bond 
distortions, site charge occupancies and spin-spin corre- 
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FIG. 8: (a) Self-consistent site charges (vertical bars) and 
spin-spin correlations (points, lines are to guide the eye) in 
the N=20 periodic ring with U = 8, V = 2.75, a 2 /K s = 1.5, 
g 2 /Kn = 0.32 for S z = 2. (b) Same parameters, except 
S z = 3. Black and white circles at the bottoms of the panels 
denote charge-rich and charge-poor sites, respectively, while 
gray circles denote sites with population nearly exactly 0.5. 
The arrows indicate locations of the defect units with nonzero 
spin densities. 



lations self-consistently with the same parameters for ex- 
cited states with S z > 0. 

In Fig. [UJa) we have plotted the deviations of the site 
charge populations from average density of 0.5 for the 
lowest S z =2 state for N = 20, for U = 8, V = 2.75 and e-p 
couplings as given in the figure caption. Deviations of the 
charge occupancies from the perfect • • • 0110- • • sequence 
identify the defect centers, as seen from Fig. [T^c) - (e). 
In the following we refer to dimer units composed of sites 
i and j as [i,j]. Based on the charge occupancies in the 
Fig, we identify units [1,2], [7,8], [13,14] and [14,15] as the 
defect centers. Furthermore, based on site populations of 
nearly exactly 0.5, we identify defects on units [1,2] and 
[7,8] as type I; the populations on units [13,14] and [14,15] 
identify them type II. Type I defects appear to be free 
and soliton like, while type II defects appear to be bound 
into a triplet state, but both could be finite size effects. 

Fig. [7|(d) and (e) suggest that the spin-spin z- 
component correlations, (SfSJ), are large and positive 
only between pairs of sites which belong to the defect cen- 
ters, as all other spins are singlet-coupled. For our char- 
acterization of units [1,2], [7,8], [13,14] and [14,15] to be 
correct therefore, (SfSj) must be large and positive if 
sites i and j both belong to this set of sites, and small 
(close to zero) when either i or j does not belong to 
this set (as all sites that do not belong to the set are 
singlet-bonded). We have superimposed the calculated 
z-component spin-spin correlations between site 2 and 
all sites j = 1 - 20. The spin-spin correlations are in 



complete agreement with our characterization of units 
[1,2], [7,8], [13,14] and [14,15] as defect centers with free 
spins. The singlet 1-1 bonds between nearest neighbor 
charge-rich sites in Fig. [Tfc) - (e) require that spin-spin 
couplings between such pairs of sites are large and neg- 
ative, while spin-spin correlations between one member 
of the pair any other site is small. We have verified such 
strong spin-singlet bonds between sites 4 and 5, 10 and 
11, and 18 and 19, respectively (not shown). Thus mul- 
tiple considerations lead to the identification of the same 
sites as defect centers, and to their characterization as 
types I and II. 

We report similar calculations in Fig.[5]Jb) for S z = 3. 
Based on charge occupancies, four out of the six defect 
units with unpaired spins in the S z = 3 state are type II; 
these occupy units [3,4], [9,10], [13,14] and [19,20]'. Type 
I defects occur on units [1,2] and [11,12]. As indicated in 
the figure, spin-spin correlations are again large between 
site 2 and all other sites that belong to this set, while they 
are again close to zero when the second site is not a defect 
site. As in the previous case, we have verified singlet spin 
couplings between nearest neighbor pairs of charge-rich 
sites. There occur then exact correspondences between 
charge densities and spin-spin correlations, exactly as for 
S z = 2. 

The susceptibility calculations in Fig. [4] and Fig. [5] are 
consistent with the microscopic calculations of spin de- 
fects presented above. As 4kp defects are added to the 
2kp background, the 2kp susceptibility peak is expected 
to broaden and shift towards higher q. This is exactly 
what is seen in the charge susceptibility, Fig. Q|a) and 
Fig. EI a) as T is increased. A similar broadening and 
shift is seen in the bond order susceptibility as well. 



V. DISCUSSIONS AND CONCLUSIONS 

In summary, the SP state in the i-filled band CTS is 
unique for a wide range of realistic Coulomb interactions. 
Even when the 4kp state is the Wigner crystal, the SP 
state can be the • • ■ 1100- ■ • BOW. For [7 = 8, for exam- 
ple, the transition found here will occur for 2 < V < 3. 
This novel T-dependent transition from the Wigner crys- 
tal to the BCDW is a consequence of the spin-dependence 
of V c . Only for V > 3 here can the SP phase be the 
tetramerized Wigner crystal 1 = = 1 • • • • • • 1 (note, 
however, that V < 4 for U — 8). We have ignored the 
intrinsic dimcrization along the ID stacks in our reported 
results, but this increases V C (U, S = 0) even further, and 
makes the Wigner crystal that much more unlikely 24 . 
Although even larger U (U = 10, for example) reduces 
V C (U,S = 0), we believe that the Coulomb interactions 
in the (TMTTF) 2 X lie in the intermediate range. 

A Wigner crystal to BCDW transition would explain 
most of the experimental surprises discussed in Section 
II. The discovery of the charge redistribution upon en- 
tering the SP phas o 15 ' 16 in X = AsF 6 and PF 6 is prob- 
ably the most dramatic illustration of this. Had the SP 
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state maintained the same charge modulation pattern 
as the Wigner crystal that exist above T 2fej? , the dif- 
ference in charge densities between the charge-rich and 
the charge-poor sites would have changed very slightly 25 . 
The dominant effect of the SP transition leading to 
1 = = 1 • • • • • • 1 is only on the charge-poor sites, 
which are now inequivalent (note that the charge-rich 
sites remain equivalent). The difference in charge den- 
sity of the charge-rich sites and the average of the charge 
density of the charge-poor sites thus remains the same^. 
In contrast, the difference in charge densities between 
the charge-rich and the charge-poor sites in the BCDW 
is considerably smaller than in the Wigner crystal 25 , 
and we believe that the experimentally observed smaller 
charge density difference in the SP phase simply reflects 
its BCDW character. 

The experiments of references [l5lfl6l should not be 
taken in isolation: we ascribe the competition between 
the CO and the SP states in X = PF 6 and AsF 6 , as 
reflected in the different pressure-dependences of these 
states^, to their having different site charge occupan- 
cies. The observation that the charge density difference 
in X = SbFg decreases considerably upon entering the SP 
phase from the AFM1 phased can also be understood if 
the AFM1 and the SP phases are assigned to be Wigner 
crystal and the BCDW, respectively. The correlation be- 
tween larger Tco and smaller Tsp^- is expected. Larger 
Tco implies larger effective V, which would lower T$p- 
This is confirmed from comparing Figs. 4 and 5: the 
temperature at which Xp(2k F ) begins to dominate over 
Xp(4k F ) is considerably larger in Fig. 4 (smaller V) than 
in Fig. 5 (larger V). The isotope effect, strong enhance- 
ment of the Tco (from 69 K to 90 K) with deutera- 
tion of the methyl groups in X = PFg, and concomi- 
tant decrease in Tsp&^£ are explained along the same 
line. Deuteration decreases u>h in Eq. . which has 
the same effect as increasing V. Thus from several dif- 
ferent considerations we come to the conclusion that the 
transition from the Wigner crystal to the BCDW that we 
have found here theoretically for intermediate V/\t\ does 
actually occur in (TMTTF) 2 X that undergo SP transi- 
tion. This should not be surprising. Given that the 1:2 
anionic CTS lie in the "weak" V/\t\ regime, TMTTF with 
only slightly smaller \t\ (but presumably very similar V, 
since intrastack intermolecular distances are comparable 
in the two families) lies in the "intermediate" as opposed 



to "strong" V/\t\ regime. 

Within our theory, the two different antiferromagnctic 
regions that straddle the SP phase in Fig. 1, AFM1 and 
AFM2, have different charge occupancies. The Wigner 
crystal character of the AFM1 region is revealed from the 
similar behavior of T N and T CO in (TMTTF) 2 SbF 6 un- 
der pressure^, indicating the absence of the competition 
of the type that exists between CO and SP, in agree- 
ment with our assignment. The occurrence of a Wigner 
crystal AFM1 instead of SP does not necessarily imply 
a larger V/\t\ in the SbFg. A more likely reason is that 
the interaction with the counterions is strong here, and 
this interaction together with V pins the electrons on al- 
ternate sites. (TMTSF) 2 X, and possibly (TMTTF) 2 Br, 
belong to the AFM2 region. The observation that the 
CDW and the spin-density wave in the TMTSF have the 
same periodicities 1 had led to the conclusion that the 
charge occupancy here is • • • 1100- • •— . This conclusion 
remains unchanged. Finally, we comment that the ob- 
servation of Wigner crystal CO 3 in (DI-DCNQI) 2 Ag is 
not against our theory, as the low T phase here is an- 
tiferromagnetic and not SP. We predict the SP system 
(DMe-DCNQI) 2 Ag to have the • • • 1100- • • charge order- 
ing. 

In summary, it appears that the key concept of spin- 
dependent V c within Eq. ([Taj) can resolve most of the 
mysteries associated with the temperature dependence 
of the broken symmetries in the ^-filled band CTS. One 
interesting feature of our work involves demonstration 
of spin excitations from the BCDW state that necessar- 
ily lead to local changes in site charges. Even for weak 
Coulomb interactions, the i-filled band has the BCDW 
character—. The effects of magnetic field on ^-filled band 
CDWs is of strong recent interes t 74 i 75 . We are currently 
investigating the consequences of the Zeeman interaction 
on the mixed spin-charge excitations of the BCDW. 
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